Glauber dynamics in a single-chain magnet: From theory to real systems 
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The Glauber dynamics is studied in a single-chain magnet (SCM). As predicted a single relaxation 
mode of the magnetization is found. Above 2.7 K, the thermally activated relaxation time is mainly 
governed by the effect of magnetic correlations and the energy barrier experienced by each magnetic 
unit. This result is in perfect agreement with independent thermodynamical measurements. Below 
2.7 K, a crossover towards a relaxation regime is observed that is interpreted as the manifestation 
of finite-size effects. The temperature dependences of the relaxation time and of the magnetic 
susceptibility reveal the importance of the boundary conditions. 

PACS numbers: 75.10.Pq, 75.40.Gb, 76.90.+d 



The design of new slow- relaxing magnetic nanosystems 
is a very challenging goal for both applications (as in- 
formation storage) and fundamental research. A well- 
known example of such systems is the single-molecule 
magnet (SMM) that shows slow reversal of the magneti- 
zation due to the combined effect of a high spin ground 
state and uniaxial anisotropy producing an energy bar- 
rier between spin- up and spin-down states [lj. When a 
magnetic field is initially applied to magnetize this sys- 
tem and then removed, the magnetization decays with a 
material-inherent relaxation time depending on the tem- 
perature. The corresponding relaxation time, r, follows 
an Arrhenius law at high temperatures and the activa- 
tion energy is equal to the barrier height, being roughly 
\D\S 2 , where D is the negative uniaxial anisotropy con- 
stant and S is the spin ground state of the molecule. At 
lower temperatures, t may saturate when quantum tun- 
neling through the barrier becomes relevant 0- 

Another research route of metastable magnetism has 
recently been explored with the synthesis of single-chain 
mag nets (SCMs) 0, 0, H 

In these materials, the slow 
relaxation of magnetization is not solely the consequence 
of the uniaxial anisotropy seen by each spin on the chain 
but depends also on magnetic correlations. The effect of 
the short-range order becomes more and more important 
when the temperature is reduced until a critical point 
is reached at T = K for ID systems. In fact, the re- 
laxation time is found to be exponentially enhanced at 
low temperatures in agreement with the pioneer work of 
R. J. Glauber devoted to the dynamics of the ID Ising 
model 0. Although it seems that there is a reason- 
able agreement between the experimental data and the 
Glauber's theory 0,0,0, we show in this communication 
that several other arguments should be considered to fill 
the gap between the theory and the experimental results. 
Firstly, it should be mentioned that the experimental sys- 
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FIG. 1: Schematic view of the chain structure 
showing the magnetic units (Mn-Ni-Mn trimers) in 
[Mn 2 (saltmen) 2 Ni(pao)2(py)2](C104)2 [4]. 



terns are not strictly Ising-like. In the simplest case, they 
are rather described by an anisotropic Heisenberg model: 



H = -2 J SiS i+1 + D SI. 



(1) 



where J is the ferromagnetic exchange constant between 
the spin units and D is the single-ion anisotropy. Sec- 
ondly, the relaxation time of each magnetic unit, intro- 
duced phcnomcnologically in Glauber's study, is a priori 
temperature dependent 7] and this argument should also 
be considered. This question is particularly important if 
each magnetic unit is by itself a slowly relaxing object. 
Finally, a real material is never perfect and, as the mag- 
netic correlation length becomes exponentially large at 
low temperature, even a very small number of defects 
should deeply affect the magnetic behavior of these one- 
dimensional systems 0. 

In this letter, we analyze the magnetic relaxation 
of a ID system, [Mn 2 (saltmen) 2 Ni(pao)2(py)2](C104)2 
(saltmen 2- = N,N'-(l,l,2,2-tetramethylethylene) 
bis(salicylideneiminate); pao - = pyridine- 2-aldoximate; 
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Py = pyridine), characterized recently (Fig. 1) |4|. At 
low temperature, this compound can be described as 
a chain of ferromagnetic coupled S — 3 [Mn iT ' r -Ni iT - 
Mn /7/ ] units. By comparing the relaxation time obtained 
from AC susceptibility with the one deduced from DC 
measurements (relaxation as a function of time), we 
show that a unique relaxation time is found over 10 
decades. We demonstrate a quantitative agreement 
between the theory and the experiment when the values 
of J and D obtained from independent thermodynamical 
measurements are considered. Moreover we show that 
finite-size effects should be considered to account for the 
crossover observed on both the magnetic susceptibility 
and the relaxation time. 

We first present a brief theoretical description of the 
expected relaxation in a SCM. Critical slowing down of 
the magnetization is expected in the vicinity of any mag- 
netic critical point. However, the ID case is singular as 
only short-range order can be observed at finite temper- 
atures. These correlations affect the dynamic response 
of the system and a slowing down of the relaxation is 
expected. The ID Ising model discussed by Glauber 
gives a single time process (i.e., a Debye relaxation) with 
an exponential enhancement of the relaxation time that 
reads (for a spin state S and with the notations of Eq. 1): 
r(T) = n{T) exp(8 J S 2 /k B T), where n{T) is the individ- 
ual relaxation time of each magnetic unit. In this rela- 
tion, the exponential factor is directly related to the tem- 
perature dependence of the correlation functions. At low 
temperatures, the corresponding energy gap for the cor- 
relation length or the magnetic susceptibility is 4 JS 2 @ . 
For any negative value of D, the correlation length re- 
mains exponentially enhanced at low temperatures [loj 
and we can still expect a Glauber relaxation. The value 
of the corresponding gap is the energy needed to cre- 
ate a Bloch wall. It remains the same as in the Ising 
limit as long as \D\ > 4 J/3 Moreover, each mag- 

netic unit exhibits uniaxial anisotropy which creates a 
barrier as in the SMM case. For magnetic units with a 
spin ground state S, we therefore expect n(T) to vary 
as Ti(T) — tq exp (\D\S 2 /fee T) when quantum tunneling 
can be ignored. Finally, we expect the relaxation time to 
vary as : 

r(T) = r exp((8J + \D\)S 2 /k B T) (2) 

The second point to consider is the influence of defects 
along the chain. Finite-size scaling of the Glauber model 
has been discussed 0, ■ In the case of open chains 
of size L, a crossover on the relaxation time has been 
predicted when the magnetic correlation length of the 
infinite chain, £, becomes of the order of L : tl(£) = 
r oo(C)/r(i/C) (where the fV is the finite-size scaling func- 
tion introduced in Ref. .12]). Following this model, the 
activation energy of the relaxation time should decrease 
from (8J+ \D\)S 2 to (4 J + |L>|)S 2 at the crossover. At 
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FIG. 2: Relaxation of the normalized magnetization from DC 
measurements. 



approximately the same temperature, we expect a satu- 
ration of x T m- 

The studied SCM is a heterometallic chain of Mn /7/ 
and Ni J/ metal ions (Fig. 1). The magnetic unit is a 
[Mn J// -Ni // -Mn /// ] trimer with strong Ni // -Mn /// anti- 
ferromagnetic interactions (</af / fes = -21 K) Q. At low 
temperature {k B T <C Jaf), the trimers can be described 
as effective S = 3 units connected within the chains by 
weaker ferromagnetic interactions (see Fig. 1). The 3D 
organization of these chains in the solid implies that they 
can be considered as magnetically isolated. The magnetic 
susceptibility of this compound becomes anisotropic with 
almost uniaxial symmetry below 50 K. Between 3.5 K 
and 7 K, a frequency dependent AC susceptibility has 
been systematically observed in both real and imaginary 
components . It approximately corresponds to a single 
relaxation time since a satisfactory fit can be obtained 
using a generalized Debye model with a small a value 
(« 0.05) Q. Below 3.5 K, the relaxation becomes too 
slow to be observed from AC susceptibility but DC mea- 
surements can still be performed. With this technique, 
we have been able to follow the relaxation down to 1.8 
K, as summarized in Fig. 2. The same shape of the 
normalized magnetization m(t) is obtained for all the in- 
vestigated temperatures, i.e. the data can be scaled into 
a single master curve. As shown in the inset of Fig. 3, 
a satisfactory fit can be obtained using a stretched expo- 
nential, i.e. m(t) = aexp(— (t/r) 13 ) + b. The b parameter 
accounts for a residual slower relaxation which only con- 
cerns a small part of the crystal (b is always less than 5% 
at any temperature). j3 is always close to 0.8 suggesting 
that some poly-dispersity may be present. The deduced 
relaxation time is depicted in Fig. 3. Note that a consis- 
tent value of t is obtained simply taking the time when 
the normalized magnetization has reached the value 1/e. 
Taken together, AC and DC data show that a crossover is 
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FIG. 3: Semi-log plot of the relaxation time r versus 1/T. 
The full and open dots were obtained from AC and DC mea- 
surements, respectively. The corresponding straight lines give 
the energy gaps: Ai/fes ~ 74K and A2//CB ~ 55K. Inset: re- 
laxation of the magnetization at 3 K. Solid line is the best fit 
obtained with the model defined in the text (top). Schematic 
view of the magnetic interactions for the model described in 
the text (bottom). 
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FIG. 4: Field dependence of the magnetization (normalized 
at the saturation value) when the magnetic field is applied 
perpendicular to the easy axis at 1.5 K. The intersection be- 
tween the two dotted lines gives the anisotropy field. Inset: 
Similar measurements than the main figure at three angles in 
the plan normal to the easy axis. The main figure is at 9 = 
0°. 



observed around 2.7 K. Indeed, above this temperature 
the dependence of the relaxation time follows an acti- 
vated law with an activation energy of Ai/fes ~ 74K\ 
Below 2.7 K, a departure from this simple behavior is 
observed and a smaller activation energy; A2/A:b ~ 55.KT 
is found around 2 K (see Fig. 3). 



Using independent magnetic measurements, we have 
determined the anisotropy parameter, D, by measuring 
the magnetization of a single crystal as a function of 
a magnetic field applied perpendicular to the easy axis 
(Fig. 4). We obtain a crossover between a linear depen- 
dence of the magnetization at low fields and a saturation 
at higher fields. Moreover, the data are almost the same 
for any direction perpendicular to the easy axis (inset 
of Fig. 4). This behavior is characteristic of the uniaxial 
anisotropy. An anisotropy field of 110 kOe is deduced 
from these data which gives D/k B = -2.5 K 13. The 
temperature dependence of the magnetic susceptibility 
is given in Fig. 5 (%T versus 1/T in a semi-log plot). 
Full dots were obtained from AC measurements at low 
frequency (i.e., below the characteristic frequency of the 
relaxation) with a Quantum Design SQUID magnetome- 
ter. Open dots were obtained from DC measurements 
with a micro-SQUID experiment 16] . The expected ac- 
tivated behavior is in fact observed above 6 K (Fig. 5). 
The deduced gap, 4JS 2 /Icb = 28K gives an estimation of 
the exchange energy J/k B = 0.78K [HIES- This gives 
(8J+ \D\)S 2 /k-B — 78K being in perfect agreement with 
the activation energy of the relaxation time, A, measured 
above 2.7 K (Eq. 2). 

Concerning the low temperature behavior, \T reaches 
a maximum around 5 K, while no anomaly was found 
on the relaxation time. This observation precludes the 
occurrence of a magnetic phase transition. On the other 
hand, this result can be explained considering finite-size 
effects (i.e., the presence of structural defects on the 
chains which limit the growth of the correlation length). 
In fact the magnetic susceptibility of a finite Ising chain 
can be calculated and leads to a saturation of xT at 
low temperature. To account for the weak decrease of 
\T observed below 5 K, weak antiferromagnetic inter- 
action, J', can be introduced between finite chains of 
n spins (as shown in inset of Fig. 3) 0|. The corre- 
sponding susceptibility can be also determined ji^l and 
gives an excellent agreement to the experimental data 
below 15 K, where the trimers can be considered as a 
5 = 3 units (Fig. 5). The deduced parameters are the 
interaction J'/fce ~ —A7mK and the average number 
of spins between two defects, n ss 110 corresponding to 
an average distance, L w 140nm. A more precise dis- 
cussion can be made using the field dependence of the 
magnetization shown at 4 K in the inset of Fig. 5. The 
same model (solid line, inset of Fig. 5) gives a reason- 
able fit with n = 90, in good agreement with the value 
deduced from the temperature dependence of the sus- 
ceptibility. As discussed previously, the finite-size ef- 
fect should also influence the relaxation time and the 
model of [13g leads to a crossover at a temperature where 
the yT product is maximum. In fact, Fig. 3 shows a 
crossover on the relaxation time, but at a smaller tem- 
perature, « 2.7K 7 compared to 5 K for the maximum 
of xT '. According to the analysis of the susceptibility 
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FIG. 5: Semi-log plot of yT versus 1/T (x is given in 
emu/mol). The full and open dots were obtained from AC 
and DC measurements, respectively. Inset: Field dependence 
of the magnetization (normalized at the saturation value) at 
4 K and with a field sweep rate of 0.14 kOe/s. Solid lines 
correspond to the fit using the expressions given in Ref. [19] 



measurements, we have generalized the model developed 
by Luscombe et al. taking into account the weak J' 
interactions [l^. The crossover on the relaxation time 
now occurs at k B T* = (4J - 4J')5 2 / ln(2n), i.e. is 
shifted at lower temperature for an ferromagnetic inter- 
action, while the activation energy below T* becomes 
A 2 /fc B = {AJ+AJ'+\D\)S 2 . To obtain T* = 2.7K, J'/k B 
has to be equal to 0.37K. This value is not in agreement 
with J'/kB — — 0.047-ftT obtained from the susceptibility. 
Nevertheless, the activation energy A2/fce ~ found 
around 2 K is close to (4 J + \D\)S 2 and supports a small 
J' value. 

In summary, the activation energy of a SCM is a com- 
bined effect of the individual kinetics of a trimer com- 
posing the chain and of the magnetic correlations. The 
energy parameters J and D deduced from the relax- 
ation time are in excellent agreement with the values in- 
dependently obtained from magnetic measurements. A 
crossover is observed at low temperature when the mag- 
netic correlation length becomes of the order of the dis- 
tance between two defects. A model of these defects in 
terms of finite-size chains coupled with weaker exchange 
interactions gives a coherent description of the temper- 
ature dependence of the magnetic susceptibility. It also 
qualitatively explains the temperature dependence of the 
relaxation time although additional ingredients like poly- 
dispersity on the chain length or exchange interaction 
may be relevant to improve the discussion of the dynamic 
properties. 
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trimer unit. It corresponds to the susceptibility of a 
chain of weakly coupled (by J ) saturated segments of 
n spins. So, in the same domain of temperatures, the 
magnetization normalized at its saturation value reads: 



m n = smh(ng[iBSH/kBT)/(sinh 2 (ng(iBSH/kBT) + 
exp^AK')) 1 / 2 . This expression has been used to fit the 
data inset Fig. 5 (solid line) . 



